Recursion method for the quasiparticle structure of a single vortex with induced 

magnetic order 
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We use a real-space recursion method to calculate the local density of states (LDOS) within a 
model that contains both d-wave superconducting and antiferromagnetic order. We focus on the 
LDOS in the superconducting phase near single vortices with either normal or antiferromagnetic 
cores. Furthermore, we study the low-energy quasiparticle structure when magnetic vortices operate 
as pinning centers for surrounding unidirectional spin density waves (stripes). We calculate the 
Fourier transformed LDOS and show how the energy dependence of relevant Fourier components 
can be used to determine the nature of the magnetic field-induced order, and predict field-induced 
LDOS features that can be tested by future scanning tunneling microscopy (STM) experiments. 

PACS numbers: 74.20.-z, 74.72.-h, 74.25. Jb, 74.25.Ha 



INTRODUCTION 



It is becoming evident that competing phases cause 
many of the anomalous properties of doped Mott insu- 
lators. An example is given by the vortex state of un- 
derdoped high-T c superconductors where antiferromag- 
netism (AF) 'pops up' near the vorticespJQ- Initial ex- 
perimental evidence for this claim came from STM exper- 
iments on YBa 2 Cu 3 O y (YBCO) and Bi 2 Sr 2 CaCu208 +I 
(BSCCO) observingweak low-energy quasiparticle peaks 
around 5-7 meV0,Q. This strongly contradicts the ex- 
pected LDOS in the vortex center of a pure BCS d-wave 
superconductor (dSC) which is dominated by the so- 
called zero-energy state (ZES), a single broad resonance 
centered at the Fermi level [5j. Further evidence for AF 
cores has come from both nuclear magnetic resonance 
measurements :Q| and muon spin rotation experiments|7J. 
The field-induced magnetization is not necessarily re- 
stricted to the core regions as determined by the co- 
herence length £. For instance, elastic neutron scat- 
tering on underdoped La2- x Sr x Cu02 (LSCO) showed 
that the intensity of the incommensurate peaks in the 
superconducting phase is considerably increased when 
a magnetic field is applied perpendicular to the Cu02 
planes or when Zn is doped into the samples 9J. Sim- 
ilar results have been found in the oxygen doped sam- 
ple La2Cu04+ a 10]. The momentum position and field- 
enhanced sharpening of this elastic signal corresponds 
to a spin density wave period of roughly eight lattice 
constants 8a extending far outside the vortex cores, sug- 
gesting that the magnetic cores operate as pinning cen- 
ters for surrounding spin density waves [ill 1 12| . This un- 
usual behavior agrees with in-field STM measurements 
on optimally doped BSCCO which found local field- 
induced checkerboard LDOS patterns with a period close 
to 4a Similar structure has been reported in zero 
field STM experiments 0| ■ Pronounced checkerboard or- 



dering has also been detected in Na-rCa9_-rCuQ9Ci9|l5j. 
More recently, Levy et al. 0] confirmed the results of Ref. 
IHi and found that the checkerboard modulation does not 
disperse with energy, and mapped out the energy depen- 
dence of the amplitude of the Fourier component corre- 
sponding to the ordering vector of the modulation. 

Theoretically, several groups have proposed that the 
origin of the unexpected behavior inside the cores is re- 
lated to locally nucleated A f1i71 fisj ] , but other scenarios 
have also been proposed^, l20j |. From a computational 
point of view, in order to model the existence of nano- 
scale inhomogeneity, it is necessary to use methods that 
easily allows one to obtain the LDOS as a function of 
energy and large real-space regions. Traditionally this 
is done by numerical diagonalization of the Bogoliubov- 
de Gennes (BdG) equations, which, at present, is typi- 
cally restricted to quite small lattices (< 40 x 40 sites). 
In this paper we use a recursion method generalized to 
the d-wave superconducting state to calculate the LDOS 
near an increasingly complex single vortex. This method 
is easily applied to large systems allowing for e.g. high- 
resolution Fourier LDOS images. First we study the pure 
dSC vortex for realistic band structure parameters rele- 
vant for overdoped cuprates. Second, we discuss the case 
of an AF vortex core in the optimally doped regime and 
focus on the spatial dependence of the expected LDOS. 
Finally, we calculate the LDOS when the vortex pins sur- 
rounding incommensurate stripe order as may be relevant 
for LSCO and underdoped BSCCO, and discuss the en- 
ergy dependence of the resulting Fourier transform. As 
opposed to most earlier theoretical work on the AF vor- 
tex problem[l7j, we focus on the final LDOS structure 
and the Fourier transformed LDOS maps which can be 
used as an STM tool to determine the nature of the field- 
induced order and the origin of the ZES splitting. Lastly, 
we compute the LDOS resulting from the recently pro- 
posed pair-density wave ordered state. 
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MODEL AND METHOD 

In the following we study the mean-field Hamiltonian 
defined on a 2D lattice 



(ij)a- i<y 

E ( a «4t4i + H - c -) + E m * (4 



(i) 



fCj-f C i{ C ii ) i 



where cj CT creates an electron with spin cr at site i, ty 
is the hopping integral to nearest (t) or next-nearest 
(t') neighbors and p is the chemical potential. The AF 
and dSC order parameters are given by njj and Ay, re- 
spectively. The Hamiltonian is the effective mean- 
field model obtained after performing two Hubbard- 
Stratonovich transformations of the extended Hubbard 
model with the on-site repulsion causing the AF, and the 
attractive nearest neighbor interaction resulting in the 
dSC. It has been used extensively in the past few years 
to gain insight into the electronic structure of p hases of 

coexisting AF and dSC order [H El HHHlUg. 

Below, we solve the Hamiltonian Q using appropriate 
Ansatze for both Ay and mj. The lack of self-consistency 
can sometimes be useful in clarifying, for instance, the 
nature/origin of vortex core states |3(|. We restrict the 
discussion to the case when the applied magnetic field 
is much smaller than H c i and consequently ignore the 
vector potential A. 

In order to obtain the LDOS near single vortices we 
use a recursion method 0i HH generalized to the super- 
conducting state. The starting point is to generate a new 
orthonormal basis of states from the recursion relation 



H \n) = a n \n) + b n+ i \n + 1) + b n \n - 1) 



(2) 



For each recursion the Greens function of the n'th level 
is generated recursively from the Lanczos coefficients a n 
and b n 



1 



h 2 r 



(3) 



Hence, the local Greens function can be found if Gn(lu) = 
for some number N, or if an appropriate analytical 
solution of Gat for an infinite chain can be attached. 
The retarded Greens function is 
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where r\ is used as an artificial smearing factor with r\ = 
0.02i. In general, it is not necessary to perform four 
recursions (I from c|^.|0); II from Cix|0); III from cJjO); 
IV from Cjf |0}) to calculate the spin-summed LDOS since 
Gil (-u) = Gi =0 (Lu) andGSH = G*l (-w). Thus, 



it is sufficient to perform only two recursions to obtain the 
total LDOS, p(i,ui), in the form of a continued fraction 

_i 

p(i, to) = Im 



ijj — <2q + irj 



u)—a\+ir} — — J-* 



+ (c^-c,a / ^ a /7 ,(6 / ) 2 ^(6 // ) 2 ), (5) 

which can be compared to the differential tunneling con- 
ductance as measured by e.g. an STM tip. Of course, 
when there is spin degeneracy (here: = 0) only one 
recursion is needed to produce the total LDOS. 

In the cases studied, we find that the Lanczos coef- 
ficients converge nicely when increasing the number of 
recursions, i.e. the system size. Below we simply per- 
form the truncation Gn — where N is some number 
of order 10 3 , and have checked that this choice does not 
affect the reported results. 



RESULTS 

In this section we use the recursion method to study 
p(i, u>) around a single vortex in the dSC state both with 
and without antifcrromagnetism in the core region. This 
is supposed to model the vortex LDOS in the overdoped 
and optimally doped regime, respectively. The core cen- 
ter is positioned at the origin (0, 0) and lengths are mea- 
sured in units of the lattice constant a. 



A single vortex without induced stripe order 

As is well-known, in an s-wave BCS superconductor 
the vortex generates states localized transverse to the 
flux line. These have been studied in great detail both 
theoretically HI |H i|| and experimentally 0. The 
core states result from the opposite sign of the super- 
current term in the particle and hole part of the BdG 
equations HHHHll. The reduction of the pair poten- 
tial near the vortex core causes only minor quantitative 
changes to these states. We have verified that the recur- 
sion method described above applied to s-wave supercon- 
ductors successfully reproduce these Caroli-de Gennes- 
Matricon bound states. 

To model an isolated d-wave vortex the following pair- 
ing potential is used 



Ay = Atanh(|r|/£)exp(iy>y), 



(6) 



where A is positive (negative) on x (y) links, r = (r^ + 
Yj)/2 and exp(i<^y) == (x + iy)/r with r = (x, y). In 
agreement with Ref. l3(M3ll we find that the suppres- 
sion of the gap in the core region results in only minor 
quantitative changes: in general the suppression tends to 
push the states slightly further toward the Fermi level. In 
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FIG. 1: LDOS along the anti-nodal (a-b) and nodal (c-d) 
direction for a single dSC vortex, (e-f) Spatial 2D structure 
of the ZES. Left column: p = *' = 0, A = 0.1*. Right column: 
fj, = -1.18*, *' = -0.4*, and A = 0.1*. 



the pure dSC state, the vortex is dominated by the well- 
known ZES@. However, the ZES is centered exactly at 
zero energy only for // = *' = 0. As opposed to the 
Caroli-de Gennes-Matricon states in the s-wave vortex, 
the ZES is made up of several states that merge to form 
the broad peak as the system size is increased in agree- 
ment with the extended nature of this peakjifj. In Fig. 
[I]we show p{i,uj) of a dSC vortex along the anti-nodal 
(a-b) and nodal (c-d) directions for A = 0.1*, £ = 5 and 
fi = t' = (a,c,e), and /x = -1.18*, *' = -0.4* (b,d,f). 
The latter parameter set provides a reasonable fit to the 
Fermi surface of slightly overdoped BSCCO with a van 
Hove singularity at u v h = —\] (4*' — /i) 2 + (4A) 2 . It is 



the (i-wave symmetry that causes the angular dependence 
(compare e.g. Fig. [IJa) andQJc)) of the higher energy 
core states at r ^ 0. The energy and amplitude of these 
core states are seen to be sensitive to the band parame- 
ters. This is also true for the ZES state as seen from Fig. 
nie-f ) : at *' = the low-energy spatial form of the LDOS 
has a star-shape due to the nodal dSC phase [il I32I |33|. 
However, for the more realistic BSCCO band parameters, 
the star is rotated with small maximum intensity along 
the anti-nodal directions, which is our prediction for the 
overdoped regime of BSCCO where competing AF order 
is expected to be absent. For *' = 0, it is well-known 
that a similar 7r/4 rotation takes place at higher ener- 
gies reve aling the spatial form of the higher energy core 
states HHIJ. 

We turn now to the simplest AF core situation where 
the suppression of the dSC gap inside the core causes a 
concomitant increase of the competing AF order 0- For 
simplicity, we model the AF core by 



m, = m(-l)( x+ ^(l - tanh(|ri|/0), 



(7) 



where = (x,y). Due to an associated local increase 
of the electron density, such vortices will in general be 
charged |3Jjj, and have been shown to remain stable when 
including the long-range Coulomb repulsion 0. Below, 
we therefore use fi = t' = 0.0 in order to model the 
close to half-filled vortex core regions as found in the 
self-consistent studies [tH. fli| . In Fig. [21 we show the 
final LDOS as a function of energy and distance to the 
AF vortex core. The induced magnetization leads to a 
splitting of the ZES as found previously fTR H^|. With 
increased magnetic order m, the resonant core states are 
pushed to higher energies and lose spectral weight. The 
vortex region is fully void of apparent core states for m > 
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FIG. 2: LDOS along the anti-nodal direction for a dSC vortex 
with antiferromagnetic core. Parameters used: fj, = *' = 0, 
A = 0.1*, m = 0.2* (a) and m = 0.5* (b). 
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FIG. 3: (Color online) DOS averaged over a coherence length 
£ within the core (solid black line), and just outside the 
core (dashed black line). The red dash-dotted line (blue 
thick dash-dotted line) show the spin-up (spin-down) resolved 
LDOS at the core center similar to the top scan in Fig. |2fa) . 



t. In this limit the low-energy LDOS has an apparent 
similarity with that of the pseudo-gap. 

Spatial averaging may mask the observability of the 
dispersive core states in Fig. For example, in Fig. |3 
we show p(i, u) (same parameters as in Fig. |2Ja)) at (0, 0) 
and (0, 7.5) averaged over a coherence length £. As seen, 
the resulting LDOS appears to be that of approximately 
non- dispersing resonant states which rapidly lose weight 
when moving away from the core region. This is similar 
to the measured differential tunneling conductance near 
the vortex cores of YBCO and BSCCO|l||. However, an 
unambiguous experimental determination of the type of 
order that induces the splitting of the ZES is important, 
and is related to the general discussion of time-reversal 
symmetry breaking for zero-energy Andreev states in d- 
wave superconductors j3||. Of course, vortices that sup- 
port AF cores would lead to a field-induced commensu- 
rate (n, 7r) signal in neutron scattering. However, purely 
from a tunneling point of view, it is possible to distin- 
guish AF and e.g. d x 2^ fU 2 + id xy induced order by us- 
ing spin polarized STM 37] . This is also shown in Fig. 
[3] where the positive (negative) bias peak is seen to be 
related to the spin-down (spin-up) LDOS, respectively. 
Importantly, this bias asymmetry will alternate with site, 
allowing for an unambiguous experimental test of AF or- 
der by a magnetic STM tip scanned through the vortex 
core region. 



A single vortex with induced stripe order 

In this section we calculate p(i,u>) around a dSC vor- 
tex which operates as a pinning center for unidirectional 
spin- and charge density modulations (stripes), expected 



to be relevant for STM experiments in the underdoped 
regime. Such inhomogeneous stripe solutions indeed ex- 
ist in a regime of intermediate AF coupling within self- 
consistent mean-field models that include the competi- 
tion between AF and dSC order^. We also briefly 
discuss the expected LDOS resulting from the recently 
proposed pair-density-wave (PDW) induced order con- 
sisting of a density wave of Cooper pairs without global 
phase coherence [sll l3f| . 

Whereas the previous section dealt with the details of 
p(i,u>) inside the core, the field-induced periodic order 
can most conveniently be studied in Fourier space 



(8) 



For site-centered stripes with spin (charge) period of 8 
(4) lattice constants, we show a typical checkerboard re- 
sult for p(i,u>) in Fig. ^a). The checkerboard pattern 
arises from including both vertical and horizontal stripes, 
which is a simple way to include the assumed slow fluc- 
tuation of the stripe domains. In addition, as shown 
recently, quenched disorder can severely smear any clear 
distinction between stripe and checkerboard symmetry 
breaking^IJ- Fig Efb) shows |p q | as a function of q x 
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FIG. 4: (Color online) (a) LDOS summed in the window lo 6 
[-0.04f, 0.04£] with m/t = 1.0. White (black) corresponds 
to high (low) LDOS. (b) jp q (tj)| at q y = 0.0 vs q x for the 
energies: u)/t = 0.0, 0.1, 0.2, 0.3, 0.4, 0.5. The curves are offset 
for clarity, (c) pq(<^) vs u) for a p-wave (A-wave, PDW) with 
t' = p = 0.0, A = 0.05£ (dashed, black lines), and the full 
vortex induced stripe situation with t' = — 0.4i, p = —1.18, 
m/t — 0.4 (red, solid line), and m/t — 1.0 (blue, dash-dotted 
line). Note that for all the curves p'q(0) = as expected for 
a ri-wave superconductor at T = 0. (d) LDOS for the clean 
dSC (PDW) shown by the dashed (solid) line. 
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for q y = 0.0 for various energies u>. The non-dispersive 
peak at the charge ordering vector Q = (27r/4,0) result- 
ing from the stripes is seen to completely dominate other 
quasiparticle interference effects. It is evident from Fig. 
EJb) that |pq(w)| displays a non-monotonic dependence 
on energy. In fact, as pointed out in Ref. for the case 
of weak translational symmetry breaking, useful informa- 
tion about the induced order and the underlying quasi- 
particle structure is contained in p'q(uj), the real part of 
Pq(uj) at the ordering vector Q. In general, for weak in- 
duced order Pq(u>) will exhibit peaks near u> = uj v h due 
to the logarithmic divergence coming from the van Hove 
points at (0, ±7r) and (±tt, 0), and near energies deter- 
mined from degeneracy points E^ = E'k+Q, where Eu is 
the quasiparticle spectrum for the homogeneous dSC[41|. 
For the simple nested Fermi surface (t! = fi = 0) and in 
the case a weak unidirectional /i-wave (fx = A sin(27r/4a:)) 
or a A-wave (A = Ao + Asin(27r/4x)), this is illustrated 
in Fig. 21c) by the black dashed lines. For the A-wave 
(/i-wave) Pq{oj) is symmetric (anti-symmetric) with char- 
acteristic peaks inside (outside) the bulk gap as well as 
weight at u> = OAt which is the van Hove energy for this 
band structure[3^ lilf . In Fig. ^c) we also show the 
full numerical result for Pq(w) in the vortex state for dif- 
ferent strengths of the magnetic order m. As seen, the 
stripe induced features in p'q(uj) are roughly symmetric 
around u> = 0.0. We find that sign changes in p'q(lo) 
at low energy to < A are only present for weak induced 
order m/t < 0.45. We have checked that Pq(lo) is de- 
termined almost entirely by the stripe order: omission of 
the vortex flow causes only minor quantitative changes 
at w < A. We expect the qualitative results presented 
in Fig. 0Jc) to apply primarily to LSCO and LBCO. In 
BSCCO, on the other hand, it is becoming clear that a 
strong component of the LDOS inhomogeneity is given 
by gap disorder E3- 

We now turn briefly to the discussion of the LDOS 
near vortices with induced PDW order which, for sim- 
plicity, is modeled with a A-wave. It is clear that PDW 
modulations cannot be the only induced order since that 
would not lead to a splitting of the ZES in the core center, 
and would not explain the enhanced spin response in the 
neutron experiments in the mixed state. Nevertheless, 
the question remains whether for certain regions of the 
phase diagram it coexists with or dominates the induced 
spin and/or charge order surrounding the cores, resulting 
in distinct features of the measured LDOS. As shown in 
Fig. 0Jc), in the case of particle-hole symmetric bands 
Pq(uj) is a good probe of the induced order since p'q(oj) 
is symmetric or antisymmetric with respect to the bias 
voltage for periodic modulations in the t\ or T3 channel 
of Nambu space, respectively. However, realistic band 
parameters and possible coexistence of other symmetry 
breakings will strongly modify /Oq(w) making detailed fit- 
ting to various assumed order parameters necessary |4l). 
Here, we propose the alternative possibility to search for 



PDW order using STM by identifying the Andreev reso- 
nant states existing in any gap modulated landscape 01- 
In Fig^Jd) we show the LDOS far away from the core 
region in the case where A is modulated by an additional 
sinusoidal wave of period four and an amplitude of 30% 
of the average gap. As seen, the Andreev states result 
in a distinct sub-gap shoulder in p(i,to) inside the bulk 
gap in regions well outside the vortex core. Such sub-gap 
structure will be approximately non-dispersing and hence 
distinct from dispersing core states extending outside the 
core region. 



CONCLUSIONS 

We have presented theoretical results for the quasipar- 
ticle structure near an increasingly complex vortex of a d- 
wave superconductor. We have discussed distinct LDOS 
features expected when magnetic or pair density wave 
order is induced by an applied magnetic field, and have 
suggested new tunneling experiments to test for field- 
induced antiferromagnetic order near the vortex cores of 
high-T c materials. 
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